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We present a determination of the shift (L) — due to the finite spatial box size L by means of A'^f = 2 
Chiral Perturbation Theory and Liischer's formula. The range of applicability of the chiral prediction is discussed. 



1. INTRODUCTION 

Finite size effects are important systematic ef- 
fects in the Monte Garlo treatment of lattice field 
theories . They become particularly large when 
the spectrum contains light particles as it hap- 
pens for QCD with light quark masses. Fortu- 
nately, chiral symmetry allows to investigate this 
region analytically by chiral perturbation theory 
(GHPT) |,|,|. In fact, through the detour of 
the chiral Lagrangian, information about infinite 
volume QGD may be obtained from finite volume 
simulations with box size L An important 
restriction is that 

F^L » 1 , M,/(4^F,) « 1 (1) 

has to be satisfied, since the chiral theory is a low- 
energy (long distance) expansion. In addition, 
one has to distinguish whether L is large com- 
pared to the Gompton wave length of the pion or 
not Q . We shall here restrict ourselves to the for- 
mer case but investigate what eq. means quan- 
titatively, by considering more than the leading 
order (LO) in the chiral expansion. 

More specifically, we study the pion mass, 
M^riL), defined as an eigenvalue of the QGD 
Hamiltonian in an L x L x L box (periodic bound- 
ary conditions), as it is extracted on a Euclidean 
lattice for sufficiently large time. The goal is 
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to calculate the shift Mt^{L) — AI-^ with Af^r 

2. LUSCHER FORMULA 

The formula that we will use , 
3 



(2) 



167r2Af^L 

dy e-V^fi+y'LF{iy) + 0(e-*^^) 



relates the finite volume mass shift to the tt-tt for- 
ward scattering amplitude F in infinite volume. 
The subleading term involves M> y^3/2 . An 
additive piece referring to the 3-particle vertex 
present in the original formula |^ does not con- 
tribute, because of parity conservation in QCD. 
Since eq. (^) builds on the unitarity of the theory, 
it holds only in full QCD, while a generalization 
to the quenched approximation seems impossible. 
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Figure 1. Integration contour in the complex 
v plane where v is the crossing variable in the 
(Minkowski space) forward scattering amplitude. 



A particular virtue of the Liischer formula is an 
economic one: inserting on the r.h.s. the one-loop 
expression in CHPT for F in infinite volume and 
integrating along the contour depicted in Fig. 1, 
one gets the asymptotic part of a formula which, 
otherwise, would have required a two-loop com- 
putation in a finite volume. 

3. INPUT FROM CHIRAL PT 

The simple a^h ^ c,d forward kinematics 
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with the crossing variable 
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means that the forward scattering amplitude, af- 
ter summing over isospin, reads 

F{iy) = A{s{v), 0, u{i^)) + 0, s{i^)) 



(5) 



where the isospin-invariant tt-tt scattering ampli- 
tude A{s,t,u) as defined in ||^ enters. 

We now use the expression for A at either lead- 
ing or next-to- leading order construct F{iy) 
via eq. (|^), and insert the result into cq. (^). At 
LO, this procedure has already been used in 



(here, M — \/2Mt^). In addition, at that order 
of CHPT, one can check explicitly that the non- 
leading terms get very small around M7rL~4. 

We now investigate, whether the subleading 
chiral corrections in A lead to a sizeable modi- 
fication of the LO mass shift eq. (0). At NLO 
the chiral expression for the scattering amplitude 
is more involved than eq. , and we evalu- 

ated the integral in eq. ^ numerically. The low- 
energy constants which are needed at NLO are 
taken from 0| , and - for the time being - the in- 
fluence of their uncertainties is not investigated. 

It is instructive to first look at the integrand 
I[y) = exp(-VM2 + y2 2,)i^(iy) in eq. (|). In 
Fig. 2 we plot it, with F constructed via LO and 
NLO CHPT, respectively. Since the calculation 
is based on CHPT and the integration variable y 
has the dimension of a mass, either the integrand 
tends to zero sufficiently fast beyond O(IOO) MeV 
or the calculation cannot be trusted. 

The two integrands are close to each other 
when the pion is light and the box size sufhciently 
large. However, in the whole range of L and Af^ 
studied in the graphs, the relative difference is 
substantial. The reason is that at NLO F{\y) 
contains contributions which grow quite strongly 
with y. To suppress these, a large L in the kine- 
matic factor exp(— Y^M^ + y^L) is needed. In- 
deed, one roughly needs F.,,L>2> and M.,r/-FV <4, 



4. MASS SHIFTS AT LO/NLO 

The LO chiral expression for the amplitude Q 

A[s,t,u)\^o^{s'Ml)lF^ (6) 



depends only on s, and henceforth the forward 
amplitude eq. is a constant and the relative 
shift 
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(7) 



4(27r)3/2 F2 (Af^L)3/2 



follows in closed form. As noticed by Gasser 
and Leutwyler, this expression agrees with their 
one-loop result |^ for Mt,:{L), when the non- 
leading terms 0(e^*^^) are dropped in the latter 
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Figure 2. The integrand I{y) at LO (dotted) and 
NLO (full line). 
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Figure 3. R={M^{L)-A'U)/A'U (with M^ = AU{L = oo)) from LO (dotted) and NLO (full line) CHPT. 



if one requires /nlo to differ only by a few percent 
from /lo- 

Note that in the graphs in Fig. 2 also values of 
Mt^L appear which are too small for the asymp- 
totic formula eq. (§) to apply. We do this in order 
to check which value of L is needed for a precise 
prediction for the integrand from CHPT. In ad- 
dition, it is an attempt to shed some light on the 
general question, what the condition Ft,L^\ in 
eq. (0) means numerically. This question is rele- 
vant also for applications of CHPT in the regime 
M^L <C 1, where eq. (^) has no basis but where 
other expansions [^J|] are applicable. 

The relative mass shift R {AU{L)- AU)/A'U 
is plotted in Fig. 3. One observes that for Ad-^L — 
4 and L = 1.5 fm the predicted mass shift is sig- 
nificantly below 1% and decreases further when L 
is increased. Even though the difference between 
LO and NLO is not so small, we expect this to 
hold beyond NLO, on the basis of our preliminary 
numerical study of NNLO effects. 

Larger mass shifts are predicted for a pion mass 
of e.g. 200 MeV, and box sizes below 2.5 fm. In 
this region, one would then like to correct lattice 
MC results for such effects, but it seems that the 
mass shifts of NLO CHPT may only be trusted 
within around 20% (of the very shift). 

We did not study L < 1.5 fm, since in this re- 
gion, the "distortion of the pion wave function" 
due to finite L may be large [|l[, and we ex- 
pect that higher order CHPT contributions will 
at most describe the onset of such effects. 

In summary, it appears that CHPT may be ap- 
plied to find out the region of parameters where 



the finite size effects are small and negligible com- 
pared to other errors of lattice QCD. We do, how- 
ever, find that in the region where a correction for 
finite size effects is necessary, a LO CHPT result 
is not sufficient to apply this correction with con- 
fidence. In order to really assess the precision of 
the NLO prediction, it will be interesting to in- 
clude yet one more order in the scattering ampli- 
tude |[7| , but then the uncertainties in the needed 
low-energy constants have to be investigated in 
detail. Such work is in progress. 
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